Breakdown of QCD factorization at large Feynman x 
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Recent measurements by the BRAHMS collaboration of high-pr hadron production at forward 
rapidities at RHIC found the relative production rate (d — Au)/(p — p) to be suppressed, rather 
than enhanced. Examining other known reactions (forward production of light hadrons, the Drell- 
Yan process, heavy flavor production, etc.), one notes that all of these display a similar property, 
namely, their cross sections in nuclei are suppressed at large xf- Since this is the region where 
X2 is minimal, it is tempting to interpret this as a manifestation of coherence, or of a color glass 
condensate, whereas it is actually a simple consequence of energy conservation and takes place 
even at low energies. We demonstrate that in all these reactions there is a common suppression 
mechanism that can be viewed, alternatively, as a consequence of a reduced survival probability for 
large rapidity gap processes in nuclei, Sudakov suppression, an enhanced resolution of higher Fock 
states by nuclei, or an effective energy loss that rises linearly with energy. Our calculations agree 
with data. 

PACS numbers: 24.85.+p, 12.40.Gg, 25.40.Ve, 25.80.Ls 



I. INTRODUCTION 

It is believed that in order to reach the smallest values 
of the light-front momentum fraction variable x in nu- 
clei, and thus access the strongest coherence effects such 
as those associated with shadowing or the color glass con- 
densate, one should go to forward rapidities, i.e., to the 
beam fragmentation region at large Feynman xp. This 
probably was the idea behind the measurements of high- 
Pt hadron production in deuteron-gold collisions at large 
rapidities performed recently by the BRAHMS collabo- 
ration [1] at RHIC. 

This proposition is based upon the usual leading or- 
der relation between x\ and xi for two colliding partons, 
x\ X2 = s/s, where s and s are the square of the cm. en- 
ergies of the colliding partons and hadrons, respectively. 
It is demonstrated in Sect. II, however, that although 
formally X2 reaches its minimal values as ii — * 1, the 
coherence phenomena vanish in this limit [2] . 

Moreover, it is shown in Sect. Ill that another effect 
causing considerable nuclear suppression for any reaction 
at large xf can be easily misinterpreted as coherence. 
The source of this suppression may be understood in 
terms of the Sudakov form factor, which is the probability 
that no particles be produced as xf — *■ 1, as demanded by 
energy conservation. Clearly, multiple interactions make 
it more difficult to satisfy this condition and therefore 
should cause an even greater suppression. 

Spectator partons, both soft and hard, may also inter- 
act while propagating through the nucleus and populate 
the large rapidity gap (LRG) that forms when xp — * 1- 
As long as production within the LRG is forbidden, the 
presence of spectators further reduces the survival prob- 
ability of LRG processes [3] . 

This mechanism of suppression can also be understood 



in terms of the Fock-state decomposition of the nucleus. 
The dominant Fock components are determined by the 
resolution of the interaction, and a nucleus can resolve 
more Fock states than a proton since the saturation scale 
Q s rises with A. Thus, one can see that the leading par- 
ton distribution involves higher multiparton Fock states 
in a nucleus and must fall more steeply toward xf = 1, 
as suggested by the Blankcnbecler-Brodsky counting rule 
[4] (see also [5, 6]). 

Note that this situation is analogous to what happens 
when one of the bound nucleons in a relativistic nucleus 
is moving with a momentum higher than the average mo- 
mentum of an individual nucleon. This is the mechanism 
of production of particles with xf > 1 in nuclear colli- 
sions [6]. 

The involvement of higher Fock states means that 
gluon bremsstrahlung is more intense in the interaction 
on a nucleus than on a proton target, i.e. leads to larger 
energy loss. Because of this, the large-x^ suppression 
may be envisioned to be a consequence of induced en- 
ergy loss. Remarkably, such an induced energy loss pro- 
portional to energy results in if scaling. This is dif- 
ferent from the energy-independent mean rate of energy 
loss found in [7-9]. The latter, however, was calculated 
with no restriction on the gluon radiation spectrum, a 
procedure that is not appropriate when xp — > I [10, 11]. 

The nuclear effects under discussion violate QCD fac- 
torization. This happens because any hard reaction is a 
LRG process in the limit xp — * 1, e.g. gluon radiation, 
is forbidden by energy conservation throughout a rather 
large rapidity interval. On the other hand, factorization 
relates this process to the parton distribution functions 
measured in other reactions, for instance deep-inelastic 
scattering (DIS), which do not have such a restriction. 
The light-cone dipole description employed here does not 
involve a twist decomposition, but apparently the lack of 
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factorization does not disappear with increasing scale. 
This will be confirmed below by explicit calculations for 
different hard reactions. 

Thus, we conclude that the nuclear suppression im- 
posed by energy conservation as x p — ► 1 is a leading 
twist effect, breaking QCD factorization. This is at vari- 
ance with the conclusion of [12], made on the basis of 
observed violation of xi scaling in J/fy production off 
nuclei and assuming validity of the factorization theorem 
[13] even in the limit Xf — > 1- 

Note that factorization is broken even for reactions on 
a free proton. For instance, the Drell-Yan reaction, which 
is a LRG process as xf — * 1, cannot be expressed in 
terms of hadronic structure functions measured in deep- 
inelastic scattering. 



II. DISAPPEARANCE OF NUCLEAR 
SHADOWING AT SMALLEST x 2 

It is convenient to study shadowing and other coher- 
ence effects in the rest frame of the nucleus. In this frame, 
the parameter controlling the interference between am- 
plitudes of the hard reaction occurring on different nu- 
cleons is the longitudinal momentum transfer, ql, related 
to the coherence length, l c = l/qx. The condition for the 
appearance of shadowing in a hard reaction is the pres- 
ence of a coherence length that is long compared to the 
nuclear radius, l c >Ra- 

Clearly at large x\ <~ 1, a hard reaction is mediated by 
a projectile valence quark and looks like a hard excitation 
of the valence quark in this reference frame. For instance, 
the Drell-Yan (DY) process looks like radiation of a heavy 
photon/dilcpton by a valence quark [14]. The coherence 
length in this case is related to the mean lifetime of a 
fluctuation q — > qll and reads [2, 14], 



lc = 



2E q a(l-a) 



k\ + (1 - a)Ml + a 2 m\ 



(1) 



Here kr and a are the transverse momentum and the 
fraction of the light-cone momentum of the quark car- 
ried by the dilepton; Mj t is the effective mass of the 
dilepton; and E q = x q sj2rajq and m q are the energy 
and mass of the projectile valence quark. This simple 
kinematic formula reflects the relation between the lon- 
gitudinal momentum transfer qh = {M 2 ^ — m q )/2Eq an d 
the coherence length l c = 1/qL- Note that the fraction of 
the proton momentum x q carried by the valence quark 
in this reference frame is not equal to x\ 1 but 
Clearly, when x\ — > 1, also a — > 1, i.e. the coherence 
length Eq. (1) vanishes in this limit, and no shadowing 
is possible. 

The onset of shadowing as function of rising coherence 
length can be approximated with good accuracy as, 



R 



shad 
A/N 



l-\a eff (T)F 2 (q L ), (2) 



Here Ra/n is the normalized ratio of the reaction rates 
on nuclear and nucleon targets. FA{q c ) may be called the 
longitudinal nuclear form factor, defined as 



Fa(Ql) 



(Ta) 



d 2 b 



dz e lz/L p A (b, z) 



(3) 



where (T A ) = 1/AJ d 2 bT\(b) is the nuclear thickness 
Ta(P) averaged over impact parameter b and is evaluated 
as the integral of the nuclear density over the longitudinal 
coordinate, T A (b) — J dzpA(b,z). The effective cross 
section 



a e ff(xi,X2,s) 



((T 2 q (ar T )) 
(<j qq (ar T )) 



(4) 



was evaluated in [2]. 

With the mean value of l c given by Eq. (3), the nucleus- 
to-deuteron ratios of DY cross sections are presented for 
different nuclei in Fig. 1 [2]. 
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FIG. 1: Shadowing in DY reaction on carbon, iron and tung- 
sten as function of X2 at = 4.5 GeV and s = 1600 GeV 2 . 
Nuclear shadowing disappears both at large and small X2, 
since the coherence length, Eq. (1), vanishes in both limits. 

As expected, in accordance with Eq. (1), shadowing 
vanishes at the smallest value of X2 that can be accessed 
in DY reactions, while according to QCD factorization it 
is expected to reach maximal strength. 

Note that the shrinkage of the coherence length to- 
wards Xf = 1 does not only lead to the disappearance 
of shadowing in this kinematic limit, but also reduces 
shadowing in DY compared to DIS throughout the en- 
tire range of x 2 . 

Note that the coherence length Eq. (1) linearly rises 
with energy, therefore the interval of x\, where the co- 
herence length contracts down to the nuclear size, shrinks 
like 1/E (see Sect. VD). 
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III. SUDAKOV SUPPRESSION, OR SURVIVAL 
PROBABILITY OF LARGE RAPIDITY GAPS 

The BRAHMS experiment [1] found a substantial nu- 
clear suppression for negative hadrons produced with 
high pt at large pseudorapidity r\ = 3.2, instead of the 
usual enhancement (see Fig. 3). Since the data cover 
rather small £2 ~ 10~ 3 , it is tempting to interpret the 
suppression as either a result of saturation [15, 16] or the 
color glass condensate [17], expected in some models [18]. 

Note, however, that the data span a region of rather 
large Xp, where all known reactions, both hard and soft, 
show considerable nuclear suppression. Moreover, avail- 
able data indicate that this effect scales with xf, rather 
than with X2, as one would expect if the scaling were the 
net effect of coherence. 

Indeed, the collection [20-22] of data depicted in Fig. 2 
for the production of different species of hadrons in pA 
collisions at different energies, with the nuclear depen- 
dence parametrized as A a , exhibit universality and Xf 
scaling. 

Data from the E866 experiment at Fermilab for nu- 
clear effects in J /iff and iff' production [23], shown in 
Fig. 8, also exhibit a strong suppression that is seen to 
scale in xp compared to lower-energy data [24] and that 
appears universal when compared with iff' [23]. Recent 
measurements of nuclear effects for J /iff production in 
D - Au collisions by the PHENIX collaboration [25] at 
RHIC are consistent with xf scaling, but they exhibit a 
dramatic violation of Xi scaling when compared with the 
#866 data [23]. 

The DY reaction is also known to be considerably sup- 
pressed at large xp [2], as one can see in Fig. 6. Unfor- 
tunately, no data sufficiently accurate to test xf scaling 
are available at other energies. 

There is a feature common to all these reactions; 
namely, when the final particle is produced with xf — * 1, 
insufficient energy is left to produce anything else. As a 
class, such events are usually called large rapidity gap 
processes. Obviously, the restriction of energy conserva- 
tion may cause substantial suppression. This is analogous 
to what happens in QED when clastic electron scattering 
occurs with no brcmsstrahlung within a given resolution; 
it is described by what is known as the Sudakov form fac- 
tor. The LRG cross section is more strongly suppressed 
as the resolution improves. 

If a large-Xi? particle is produced, the rapidity interval 
to be kept empty is Ay = — ln(l— xf). We describe parti- 
cle production via perturbative gluon radiation [26] with 
subsequent nonperturbative hadronization. Assuming as 
usual an uncorrelated Poisson distribution for gluons, the 
Sudakov suppression factor, i.e. the probability to have 
a rapidity gap Ay, becomes 

S(Ay) = e -^ lG(A ^ , (5) 

where in our case na(Ay) is the mean number of gluons 
that would be radiated within Ay if energy conservation 
were not an issue. 



Note that even in the case where no gluon is radiated 
within the rapidity gap, the hadronization can easily fill 
the gap with particles. The probability that this does 
not happen is another suppression factor which, however, 
is independent of target and cancels in the nucleus-to- 
proton ratio. 

The mean number (nc(Ay)) of gluons radiated in the 
rapidity interval Ay is related to the height of the plateau 
in the gluon spectrum, (na(Ay)) = Aydnc/dy. Then, 
the Sudakov factor acquires the simple form, 

S(x F ) = (1 - x F ) dnG/dy ■ (6) 

The height of the gluon plateau was estimated by Gunion 
and Bertsch [26] as, 

The value of a 8 was fitted [26] to data on pion multiplicity 
in e + e~ annihilation, where it was found that a s = 0.45. 
This is close to the critical value a s = a c = 0.43 [27] 
and to the value (a s ) = 0.38 calculated within a model 
of small gluonic spots when averaged over the gluon ra- 
diation spectrum [28]. For further calculations we take 
a s = 0.4, which gives with high accuracy dna/dy = 1, 
i.e. the Sudakov factor, 

S{x F ) = 1 -x F ■ (8) 

Amazingly, this coincides with the suppression factor ap- 
plied to every additional Pomeron exchange in the quark- 
gluon string [29] and dual parton [30] models based on 
the Regge approach. 

Clearly, on a nuclear target, the Sudakov suppression 
factor should fall more steeply as xp — ► 1 since multi- 
ple interactions enhance the transverse kick given to the 
projectile parton and therefore tend to shake off (i.e. to 
radiate) more gluons. This can be understood in terms of 
the Fock state decomposition. Specifically according to 
the counting rules [4], the behavior of the single parton 
distribution function for x\ — > 1 depends on the number 
of constituents in the particular Fock state. A nucleus 
having a higher resolution, controlled by the saturation 
scale Q s [16, 17], resolves more constituents and thus 
results in a steeper fall off of the distribution function 
toward x\ = 1. 

We come to the nontrivial conclusion that the effective 
parton distribution function in the beam hadron depends 
on the target. Such a process-dependence constitutes 
an apparent breakdown of QCD factorization and is a 
leading twist effect. 

One can also formulate this suppression as xf — ► 1 
as a survival probability of the LRG in multiple inter- 
actions with the nucleus. Clearly, every additional in- 
elastic interaction contributes an extra suppression fac- 
tor S(xf). The probability of an n-fold inelastic colli- 
sion is related to the Glauber model coefficients via the 
Abramovsky-Gribov-Kancheli (AGK) cutting rules [19]. 
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Correspondingly, the survival probability at impact pa- 
rameter b reads, 



exp[-<7-: 

A 



:N T A (b)] 

x [o% r T A (b)S n (x F )] n . (9) 

it. 

n=l 

In this expression particles (gluons) are assume to be 
produced independently in multiple rescattering, i.e. in 
Bethe-Heitler regime. Of course it should be corrected 
for effects of coherence which turn out to be either small 
or absent, i.e. can be neglected. Indeed, at small xf the 
Sudakov factor S(xf — > 0) — ► 1, and Eq. (9) takes the 
form of the standard Glauber expression for absorptive 
hadron-nucleus cross section. In this case the coherence 
effects are known as Gribov inelastic shadowing correc- 
tions which are known to be quite small, few percent [31]. 

In another limiting case xf — ► 1 energy conservation 
allows only radiation of low-energy gluons having short 
coherence time. Therefore, particles are produced inco- 
herently in multiple interactions, and Eq. (9) is legiti- 
mate. 

At large xf ~ 1, Eq. (9) is dominated by the first term; 
therefore, integrating over impact parameter, one gets for 
the nucleus-to-proton ratio, Ra/ p (xf — > 1) ~ A 1 / 3 . This 
expectation is confirmed by a measurement [32] of the 
A-dependence of the cross section for the LRG process 
pA — > pX, quasi-free diffractive excitation of the nucleus. 
The single diffraction cross section was found to be 



V. 



pA 
diff 



1 

0.925 



dxi 



da(pA — > pX) 
dxF 



(10) 



with a = 0.34 ± 0.02, consistent with the above expecta- 
tion. 



IV. PRODUCTION OF LEADING HADRONS 
WITH SMALL p T 

The collection of data from [20-22] for the production 
of different species of particles in pA collisions, depicted 
in Fig. 2, exhibits quite strong and universal nuclear sup- 
pression at large xf- Moreover, these data spanning the 
lab energy range from 70 to 400 GeV, demonstrate that 
the nuclear effects scale in xp- 

It is natural to relate the observed suppression to the 
dynamics discussed above, which is close to the descrip- 
tion of soft inclusive reactions within the quark-gluon 
string [29], or dual parton [30] models. 

The nuclear effects can be calculated summing oven n 
and integrating over impact parameter in Eq. (9), 

R A/N (x F ) = - \ j / d 2 be~^ T ^ 

(1 - XppabsA J 
,(l-x F )cr abs T A {b) _ i _ ^-Qj 




0.2 0.4 0.6 0.8 1. 



FIG. 2: The exponent describing the ^-dependence (oc A a ) 
of the ratio for the production of different hadrons in p — Au 
relative to pp collisions as function of xf. The collection of 
data and references can be found in [20-22]. The curve is a 
result of a parameter free calculation using Eq. (11). 



In the Glauber model, the effective cross section is the 
familiar inelastic NN scattering cross section. However, 
the actual number of collisions, which determines the 
value of the effective absorption cross section, is sub- 
ject to considerable modification from Gribov's inelastic 
shadowing corrections, which make the nuclear medium 
substantially more transparent. These corrections con- 
siderably reduce both the number of collisions and the 
effective absorption cross section a a bs [31, 33]. 

In order to compare with data, the nuclear effects are 
parametrized as Ra/n oc A a , where a varies with A. We 
use A=40, for which the Gribov corrections evaluated in 
[33] lead to an effective absorption cross section a a b s ~ 
20 mb. Then the simple expression Eq. (11) explains the 
observed xf scaling and describes rather well the data 
depicted in Fig. 2. 

One may wonder why a(x F ) plotted in Fig. 2 does 
not reach values as small as 1/3, even when x F — > 1- 
As mentioned, this exponent varies with A, and simple 
geometrical considerations may be accurate only for suf- 
ficiently heavy nuclei. In the case of diffraction [33] , there 
is a specific enhancement of the effective absorption cross 
section that makes a smaller in this case, in good accord 
with data [32]. 

Note that our description is very close to that in the 
dual parton model (or quark-gluon string model) [34]. 
However, we present a different interpretation of the same 
phenomena and introduce Gribov corrections for inelas- 
tic shadowing, which substantially reduce the number of 
collisions. 
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V. HIGH-pt HADRON PRODUCTION AT 
FORWARD RAPIDITIES, THE BRAHMS DATA 

The BRAHMS collaboration performed measurements 
of nuclear effects for production of negative hadrons at 
pseudorapidity r\ = 3.2 and transverse momentum up to 
Pt ~ 4 GeV. Instead of the usual Cronin enhancement a 
suppression was found, as one can see from Fig. 3. 
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FIG. 3: Ratio of negative particle production rates in d — Au 
and pp collisions as function of pr- Data are from [1], solid 
and dashed curves correspond to calculations with the diquark 
size 0.3 fm and 0.4 fm respectively. 

First, consider the rather strong suppression of the 
data at small px- One can understand this in terms of 
the simple relation for the pT-hitegrated cross sections, 



/ 



d p T 



da I dn 

di]d 2 pT \ drj 



(12) 



The mean number of produced particles per unit rapid- 
ity {dn/drj) has an A-dependence which varies with ra- 
pidity. Particle production at mean rapidity is related 
to the radiation of gluons, whose multiplicity rises as 
dnc/drj oc A 1 ! 3 (for the moment we neglect gluon shad- 
owing and assume the Bethe-Heitler regime for gluon ra- 
diation). Since the inelastic cross section af n oc A 2 / 3 , the 
integrated inclusive cross section, Eq. (12), rises linearly 
with A. This is in accordance with the AGK cancellation 
[19] of shadowing for the inclusive cross section known as 
the Kanchcli-Muellcr theorem. 

Nuclei modify the px distribution of radiated gluons, 
an effect known as the color glass condensate (CGC) 
[16, 17] or Cronin effect [35, 36]. Due to this effect, glu- 
ons are suppressed at small px, enhanced at medium pt, 
and are unchanged at large px- Gluon shadowing, or the 
Landau-Pomeranchuk effect, is a part of the CGC and re- 
duces the total number of radiated gluons more strongly 
at small than at large px- Thus, the observed strong sup- 



pression of small pr particle production at mid rapidities 
is a manifestation of the CGC. 

One, however, should be careful with the interpreta- 
tion of data in terms of the CGC, which is supposed to 
be a result of coherence between different parts of the 
nucleus. It turns out that nuclear modifications of the 
transverse momentum distribution occur both in the co- 
herent and incoherent regimes. While the former can 
be an effect of the CGC, the latter has little to do with 
this phenomenon. In particular, the RHIC data at mid 
rapidities are in the transition region, i.e. particles are 
produced coherently on the nucleus at small px < 1 GeV, 
but incoherently at larger px [36] . 

The suppression at small-py observed at rj ~ 3.2 is 
even stronger than at mid rapidities. At this rapidity, 
the overall scale of the suppression is related to the fact 
that particle production is dominated by fragmentation 
of the projectile valence quarks. Gluons are additionally 
suppressed due to softness of the gluon fragmentation 
function leading to a substantially larger value of x\ for 
gluons than for pions. Therefore, the origin of the the 
suppression is quite different from that at mid rapidity 
[37]. Because the number of valence quarks is fixed and 
equal to three when integrated over rapidity (Gottfried 
sum rule), the number of valence quarks produced with 
xp — > 1 must be even smaller, and accordingly the ra- 
tio of the pT-integrated inclusive cross sections should be 
suppressed well below unity. In this case we can use 
either our results or the data plotted in Fig. 2, both 
of which suggest a suppression factor of approximately 
A~ - 3 w 0.2, in good agreement with the BRAHMS data. 
This suppression is not affected either by the CGC or 
gluon shadowing. 

Note that the dominance of valence quarks in the pro- 
jectile proton leads to an isospin-biased ratio. Namely, 
negative hadrons with large px close to the kinematic 
limit are produced mainly from u, rather than d, 
quarks. Therefore, more negative hadrons are produced 
by deuterons than by protons, and this causes an en- 
hancement of the ratio plotted in Fig. 3 by a factor of 
3/2 [37]. Further on, we take care of this by using proper 
fragmentation functions for negative hadrons. 

The cross section of hadron production in dA(pp) colli- 
sions is given by a convolution of the distribution function 
for the projectile valence quark with the quark scattering 
cross section and the fragmentation function, 



da 



d 2 px drj 
da[qA(p)} 



1 

J dz fq/d( v ){xuq T ) 



d 2 qx drj 



Here 



q T= p T / z 
IT „ 



Dh-/ q {z) 



(13) 



(14) 



We use the LO GRV parametrization [38] for the quark 
distribution in the nucleon. As we explained above, the 
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interaction with a nuclear target does not obey factor- 
ization, since the effective projectile quark distribution 
correlates with the target. Summed over multiple inter- 
actions, the quark distribution in the nucleon reads, 



q/N< 

fd 2 b 



fqJlV ( x U<1t)= C fq/N {Xl , <7t) 

•ia eff T A {b) _ e -ffe//T A (6)l 



(l- Xl )Jd 2 b [1 - e -<W/T A (6)] 



(15) 



Here the normalization factor C is fixed by the Gottfried 
sum rule. 

The cross section of quark scattering on the target 
in Eq. (13) is calculated in the light-cone dipole ap- 
proach [39, 40], which provides an easy way to incorpo- 
rate multiple interactions. Obviously, the pr distribution 
of hadrons in the final state is affected by the primordial 
transverse motion of the projectile quarks. In our cal- 
culation, we separate the contributions characterized by 
different initial transverse momenta and sum over three 
different mechanisms of high-py production. 



B. Diquark break up qq — > qq. 

At larger kx the interaction resolves the diquark, so 
its break-up should be considered. This contribution is 
calculated in accordance with Refs. [42, 43], 

dajqqA -> qX) _ f d 2 nd 2 r 2 4 fc T ( r -_ r - 2 ) 



d 2 k T d 2 b J 2(2tt) 2 

x *k(n) ^ D {r 2 ) [2 - e-W^O^W 

-K"(*V2)T.4(b) _ 



e -I ff ~(r 2 /2)T A (fc) 



(17) 



The diquark wave function is also assumed to have a 
Bessel function form, with a mean quark separation of 
0.2 — 0.3 fm. There is much evidence that such a small 
diquark represents the dominant quark configuration in 
the proton [41]. 



A. Quark-diquark break up of the proton. 

We employ the quark-diquark model of the proton with 
the ud diquark small compared to the proton radius [6, 
41]. Correspondingly, the third valence quark external 
to the diquark has much smaller transverse momentum 
than the two others sitting inside the diquark. As a first 
mechanism, we consider proton breakup p — > qq + q. We 
treat the diquark {qq} as point-like and integrate over its 
momentum. Then the kr distribution of the projectile 
valence quark, after propagation through the nucleus at 
impact parameter b, is given by [42, 43], 



da(NA -> qX) 
d 2 k T d 2 b 



d 2 rid 2 r 2 



3 ifc r (ri-r 2 ) 



1 + e -^g(Pi-^)T A (b) 



-*er£(ft)IU(6) 



(16) 



Here the quark-diquark wave function of the nucleon is 
taken in a form that matches the known perturbative 
QCD behavior at large transverse momenta, *jv(r) oc 
Ko(r/Rp), where Kq is the modified Bessel function, 
Rp = I r 2 h , and r c h is the mean charge radius of the pro- 
ton. We assume that the quark's longitudinal momentum 
dependence factorizes and is included in f q /N( x i)- The 
dipole cross section <r^ q {r) is taken in the saturated form 
[43] , inspired by the popular parametrization of Ref. [44] , 
but adjusted to the description of soft data. 

This contribution dominates the low transverse mo- 
mentum region kr < 1 GeV. 



C. Hard gluon radiation q — > Gq. 

At large kx, the dipole approach should recover the 
parton model [45], where high momentum transfer pro- 
cesses occur (in leading order) as binary collisions with 
the transverse momentum of each final parton of or- 
der kx- Clearly, this is different from the description 
in Eqs. (16)-(17), where one assumes that the projectile 
valence quark acquires high transverse momentum as a 
result of multiple rescatterings, while the radiated glu- 
ons that balance this momentum are summed to build 
up the dipole cross section. The latter is fitted to DIS 
data involving gluons of rather low transverse momenta. 
Therefore, one should explicitly include in the dipole de- 
scription radiation of a gluon with large transverse mo- 
mentum that approximately balances kr, i.e. the process 
qN — > qGX. In the dipole approach, the cross section is 
given by the same formula, Eq. (16), except that the nu- 
cleon wave function is replaced by the quark-gluon light- 
cone wave function, ^n{tt) ^qG( r T) [42, 46], where 
[43] 



2i [as r T ■ e* 
*^ (rr ) = -_W_ exp 



(18) 



and r = 0.3 fm. Such a small mean quark-gluon sepa- 
ration is a result of a phcnomenological analysis of data 
for soft single diffraction pp ~ * pX. The only way to 
explain the abnormally small triplc-Pomcron coupling, 
which is translated into very weak diffractivc gluon radi- 
ation, is to assume that the Weizsacker- Williams gluons 
in the proton are located within small spots [49]. The 
spot size r was fitted to diffractive data [43], and the 
result r = 0.3 fm agrees with both lattice calculations 
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[50] and also with the phenomenological model of the 
instanton liquid [51]. 

Notice that at small kx there is a risk of double count- 
ing in such a procedure, since the radiated gluon may be 
counted twice, explicitly and implicitly as a part of the 
dipolc cross section. However, the first two contributions, 
Eqs. (16)-(17), and the last one dominate in different re- 
gions of kr and we found their overlap very small. 



D. Gluon shadowing 

Although the BRAHMS data involve rather large val- 
ues of Xf, the corresponding values of x 2 are so small 
that the considerations of Sects. II and VII lead to little 
reduction of the coherence length. Nevertheless, gluon 
shadowing corrections are expected to be quite small 
even at very small x 2 (e.g. see predictions for LHC in 
[36]). This is related to the presence of small size gluonic 
spots in nucleons, discussed in the previous Sect. V C. In- 
deed, one goes to small x 2 in order to make the coherence 
length longer (l c oc l/x 2 ) and thus to arrange an overlap 
between the parton clouds which belong to different nu- 
cleons separated in the longitudinal direction. However, 
even if this condition is fulfilled, the gluon clouds hardly 
overlap in impact parameter, if they are shaped as small 
size spots. Expecting the strongest gluon shadowing at 
Q 2 — > one can estimate the shadowing correction to be, 



Ra{b) ~ 1 — exp 



~rlT A (b) 



(19) 



For heavy nuclei this estimate gives a rather weak shad- 
owing, about 10%, in good accord with more accurate 
calculations [43, 66] or with the NLO analysis of nuclear 
DIS data [47]. 

We calculate gluon shadowing corrections within the 
dipole approach with the light-cone gluon distribution 
function, Eq. (18). The details can be found in [43, 66]. 



E. Comparison with data 

First of all, one should confront the model with the 
PT-dependent cross section of hadron production in pp 
collisions. Although the nuclear effects under discussion 
are not sensitive to this dependence, which mostly cancels 
in the dA/pp ratio, this would be a stringent test of the 
model. Our calculations are compared with pp data from 
the BRAHMS experiment at 77 = 3.2 in Fig. 4. In view 
of the isospin asymmetry of leading particle production 
mentioned above, it is important to use proper fragmen- 
tation functions. The standard ones extracted from data 
on e + e~ annihilation give a sum of positive and negative 
hadrons, while one needs only fragmentation functions 
for production of negative hadrons in order to compare 
with the BRAHMS data. We use these negative frag- 
mentation functions in our calculations [48]. 
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FIG. 4: Number of negative hadrons versus pt produced in 
pp collisions at y/s = 200 GeV and pseudorapidity r\ — 3.2. 
Our calculations, given by the solid curve, are compared with 
BRAHMS data [1]. 



Now we are in a position to predict nuclear effects 
employing the dipolc formalism and the mechanism de- 
scribed above. The results are compared with the 
BRAHMS data for the minimum-bias ratio[l] in Fig. 3. 
One can see that this parameter-free calculation docs 
not leave much room for other mechanisms, including 
a strong CGC. On the other hand, our calculations do 
include the CGC via explicit gluon radiation and via 
gluon shadowing. This is, however, a rather moderate 
effect due to smallness of the gluonic spots in nucleons 
[43, 49, 51], as is described by the quark-gluon wave func- 
tion in Eq. (18). 

It is interesting also to check whether the predicted de- 
pendence of the ratio on impact parameter is supported 
by data. Our results for the ratio of the cross sections 
of central and semicentral to peripheral collisions are de- 
picted in Fig. 5 in comparison with BRAHMS data [1]. 
The agreement is rather good. 

Although the BRAHMS collaboration has presented 
their results for nuclear effects as a function of rapidity, 
we skip this comparison, since it cannot fail. Indeed, 
predictions [36] for the Cronin effect at r\ = published 
in advance of data were quite successful. That region is 
dominated by production and fragmentation of gluons. 
The very forward region, which is under consideration 
now, is dominated by production and fragmentation of 
valence quarks, and our calculations are successful here 
as well. For these reasons, we should not be much off 
the data at any other (positive) rapidity, which differs 
only in the relative contributions of valence quarks and 
gluons. 
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FIG. 5: Ratio of negative particle production in central (0- 
20%) and semi-central (30-50%) to peripheral (60-80%) d-Au 
collisions, shown by closed and open points respectively. The 
results of the corresponding calculations are depicted by solid 
and dashed curves. 



VI. NUCLEAR SUPPRESSION OF DILEPTONS 
AT LARGE x F 



A. Rest Frame Description 

According to the rest frame interpretation [14] (see also 
[42, 55-58]), the DY process looks like fragmentation of 
a projectile quark into a dilcpton via bremsstrahlung of 
a heavy photon. One can cither calculate this pertur- 
batively and use the quark distribution function at the 
corresponding scale, or employ the soft quark distribu- 
tion functions and use a phcnomcnological fragmentation 
function q —* qll. The latter approach, already used in 
[2], is appropriate here, and we will apply it using the 
cross sections of soft production of valence quarks in pp 
and pA collisions fitted in the following section. To the 
extent that these cross sections are subject to nuclear 
suppression at large xp, the DY process should be sup- 
pressed as well. 

We perform calculations within the same formalism 
used in [2], but the source of suppression is not a simple 
initial state energy loss, but an effective one that results 
from the nuclear modification of the Fock state decom- 
position discussed above. The nucleus-to-deuterium ratio 
reads, 



Ra/d{x\) 



2 j dx ^pn Dn/qif) 

Xl 



(20) 



It was first observed in the E772 experiment at Fermi- 
lab [52] that the DY process is suppressed at large Xp. 
Two mechanisms that may possibly be responsible for 
this effect have been considered so far, energy loss in the 
initial state [2, 53, 54] and shadowing [2, 52]. To make 
the interpretation more certain, and to disentangle these 
two options, one can select data with the dilepton effec- 
tive mass sufficiently large to assure that the coherence 
length is too short for shadowing. An example of such 
data for the ratio of tungsten to deuterium is depicted in 
Fig. 6. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 
*1 

FIG. 6: Ratio of DY cross sections on tungsten and deuterium 
as a function of x\, at large dilepton masses to eliminate nu- 
clear shadowing. 



Here we implicitly take into account the difference in the 
isospin composition of deuterium and tungsten. We also 
incorporate the contribution of projectile antiquarks and 
target quarks using the STEQ quark distribution func- 
tions [59]. 



B. Nuclear suppression of valence quarks 

To evaluate Eq. (20) , one needs to know the cross sec- 
tions of soft valence quark production in pp and pA colli- 
sions. To obtain this, we turn the problem around, trying 
to be more model independent, and get the nuclear sup- 
pression of valence quark jets directly from data. For this 
purpose, we fitted data [20] for pion production in pp and 
pA collisions at 100 GeV. We describe the spectrum of 
produced pions as, 



da(pA — > ttX) 



dx 



dx\ 
da(pA 



qX) 



dx 



D„ /q (x/ Xl ) . (21) 



The fragmentation function q — > ir, D 7T f q (z), is known. 
We use the form suggested by the Regge approach [29], 
D n+/U = D n - /d = (1 - z)-<**+ x and D n+/d = D n - /u = 
(1 - z)D w+/u , where A « 1/2. 

The unknown function in Eq. (21) is the cross sec- 
tion for quark jet production in pp or pA collisions. 
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We parametrized these cross sections by a simple x- 
dependence, oc (1 — x c ) u / ' \fx and oc (1 — x e ) d / ' \fx for pro- 
duction of u and d quarks respectively, performed a fit to 
Fermilab data [20] at 100 GeV, and found u = 1.85±0.07, 
d = 3.05 ± 0.15 for pp collisions; u = 2.00 ± 0.11, 
d = 4.15 ± 0.33 for p-Al collisions; u = 2.03 ± 0.14, 
d = 4.00 ± 0.33 for p-Pb collisions. The quality of the fit 
can be seen in Fig. 7. 




10 r 10 r 

0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 
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FIG. 7: Fit to data [20] for pp(A) — > it with parametrization 
described in text. The results are shown by dashed and solid 
curves for pp and and p — Pb collisions respectively. 

While the parameter u does not show any strong A- 
dependence, the value of d rises with A making the in- 
dependence of the cross section steeper. This is no sur- 
prise, since d quarks have a steeper distribution func- 
tion in the proton. Unfortunately, the data are not suf- 
ficient to fix well the parameter e, which is found to be 
e = 1.5 ±0.9. 



C. Comparison with data 

Droll Yan results based on Eq. (20) are compared with 
the E772 data in Fig. 6. For D u/q (z,M 2 ) we used the 
phenomenological fragmentation function Dj l / q (z, M 2 ) oc 
(1 — z) 3 that was fitted to data in Ref. [2] and found to 
be independent of M 2 within errors, and for the quark 
production cross sections we used the fit performed in the 
previous section. The parameter e, which is poorly de- 
fined by the data, affects the x\ dependence of DY cross 
section ratio but hardly varies the amount of suppression 
at large Xp. We select e = 1.3, which is well within the 
errors and provides for Ra/d{ x i) a shape similar to the 
data. 

Note that this description of nuclear effects appears 



different from the energy- loss scenario employed in [2]. 
There is, however, a strong overlap between the two 
mechanisms. The nuclear modification of the projectile 
Fock decomposition is just a different way of calculat- 
ing nuclear effects due to gluon bremsstrahlung, which 
is the source of the induced energy loss. Therefore the 
current approach might be called an effective energy-loss 
description. The results, however, are different. Namely, 
the simple mean energy loss of a single parton leads to 
nuclear modifications that scale in Axi, i.e. vanish with 
increasing energy at fixed x\. In contrast, the current 
multi-parton effective energy loss rises linearly with en- 
ergy and leads to an x\ scaling in good accord with data. 

Note that our interpretation of nuclear effects in the 
DY process is quite different from the description in [34] , 
where nuclear effects are predicted to scale in x 2 and 
are possible only if the coherence length is longer than 
the nuclear size. On the contrary, we expect the nuclear 
suppression to scale in xp. 



VII. CHARMONIUM SUPPRESSION AT 
LARGE x f 

Nuclear suppression of Charmonium production has 
been observed to be steeply increasing at large xf ~ 1 in 
many experiments. Understanding this effect has been 
a challenge for a long time. Although the first data [24] 
were well-explained and even predicted [60] by an energy- 
loss mechanism [53], later data on J/^ production at 
higher energies demonstrated that this mechanism is not 
sufficient, since it does not explain the observed xp scal- 
ing. 

This problem was studied within the dipole approach 
in [61]. A substantial part of the suppression was found 
to be a higher twist effect related to the large size of char- 
monia. Such a suppression is frequently identified with 
simple final-state absorption. However, at high energies 
a cc fluctuation of a projectile gluon propagates and at- 
tenuates through the entire nucleus [62-64]; moreover, 
in the case of hadroproduction, such a dipole is colored. 
A description of this process in terms of the light-cone 
dipole approach was developed in [61], and the cross sec- 
tion was calculated in a parameter-free way. 

The rest of the suppression observed experimentally 
was prescribed to be the effect of energy loss and gluon 
shadowing. Within this interpretation of data, the ob- 
served x\ scaling looked like an accidental compensa- 
tion of energy-loss corrections decreasing with energy and 
gluon shadowing effects rising with energy. Correspond- 
ingly, an approximate x 2 scaling (broken at low energies 
by energy loss) was predicted in [61] for energies ranging 
between Fermilab (fixed targets) and RHIC. However, 
the recent data from the PHENIX experiment at RHIC 
found a dramatic violation of x 2 scaling in strict contra- 
diction with this prediction. In fact, any parton model 
based on QCD factorization predicts x 2 scaling and con- 
tradicts this data. 
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We think that the higher twist nuclear shadowing was 
correctly calculated in [61], but that the gluon shadow- 
ing was miscalculated and led to the incorrect predictions 
for RHIC [61, 65]. What was missed in [61] is the shrink- 
age of the coherence length towards the kinematic limit. 
This effect, found for the DY reaction in [2] and discussed 
above in Sect. II, is even more important for gluon shad- 
owing. Due to proximity of the kinematic limit for J/\f' 
production with xp — > 1, the effective value of x 2 is sub- 
stantially increased, 



X2 



x 2 



1 — X\ 



(22) 



Additionally, the coherence length available for gluon 
shadowing gets another small factor P G , 



: pG 



pG 



wln x 2 M^rriN 



Xi(l - Xi) 



(23) 



This factor was evaluated in [66] as P G « 0.1. Thus, 
we conclude that for the kinematics of the E772/E866 
experiments the coherence length for gluon shadowing 
does not exceed 1° < 0.8 fm, i.e. no gluon shadowing is 
possible. Therefore, one should search for an alternative 
explanation of the data. 

Obviously, the same Sudakov effect that causes the 
large xp suppression of other particles, in particular light 
hadrons and lepton pairs, affects the charmonium pro- 
duction as well. We can use the same results for nuclear 
softening of the produced valence quark, fitted to data, 
as for the DY reaction. Nevertheless, the phenomenolog- 
ical fragmentation function q — > tyq should be different. 
Fitting data for J/^ production in pp collisions we found 
D<n/ q {z) oc (1 — z) 1 - 6 . Then, using the same convolution 
of the distribution function of the produced quark and 
the fragmentation function as in (20), we arrive at the 
suppression depicted by the dashed curve in Fig. 8. 

After the higher twist shadowing calculated in [61] is 
added, the result is closer to the data, as shown by the 
solid curve in Fig. 8. We see that the scale of the nu- 
clear suppression at large xp agrees with the data, al- 
though the shape of the xp dependence needs to be im- 
proved. This problem is a consequence of the oversimpli- 
fied parametrization for the distribution function of the 
produced quark. Unfortunately, the errors of the avail- 
able data for light hadron production are too large and 
do not allow use of more sophisticated parametrizations. 

Since the Sudakov suppression scales in xi, one should 
expect an approximate xp scaling, which is indeed ob- 
served in data. Therefore a similar xp dependence is 
expected at RHIC and LHC, but the onset of x 2 scaling, 
which has been naively expected at high energies [61], 
will never occur. Preliminary data from RHIC [25] con- 
firm this. Additionally, they do not show any apprecia- 
ble effect of gluon shadowing, in spite of smallness of x 2 ■ 
This agrees with the finding of [43] , that the leading twist 
gluon shadowing is very weak. There is the possibility of 
a strong higher twist effect that may make gluon shad- 
owing in the charmonium channel stronger than in DIS 




FIG. 8: Tungsten-to-beryllium cross section ratio for charmo- 
nium production as function of x\. Data are from the E866 
experiment [23] . The dashed curve shows the contribution of 
the extra Sudakov suppression extracted from data for soft 
hadron production. The solid curve also includes the higher 
twist shadowing related to the nonzero cc separation. 



[61]. Such a prediction is based, however, on an ad hoc 
phenomenological potential model and may be incorrect. 

Open charm production is expected to have a simi- 
lar Sudakov nuclear suppression, like that shown by the 
dashed curve in Fig. 8. The higher twist shadowing re- 
lated to the nonzero cc separation does exist [67], but 
is much weaker [67]. The leading twist gluon shadowing 
is also rather weak even at the RHIC energies. As for 
higher twist corrections, in the potential model [67] they 
may be large, but as we mentioned, this is not a solid 
theoretical prediction. 

Note that our description of nuclear suppression of 
heavy quarkonia is quite different from the model pro- 
posed in [34]. That model involves three unknown pa- 
rameters fitted to the nuclear data to be explained. The 
key parameter is the absorption cross section for a dipolc 
consisting of a colored heavy quark pair QQ and light 
quarks. However, the pair of heavy quarks that eventu- 
ally forms the detected quarkonium, and the comoving 
light quarks, cannot "talk to each other" due to Lorentz 
time dilation during propagation through the nucleus. 
In other words, it does not make any difference whether 
the accompanying light quarks are primordial or are cre- 
ated during hadronization of the color-octet QQ pair. 
Therefore, the multiple interactions of such a large dipole 
should not be treated as absorption for production of a 
colorless QQ pair via color neutralization. In contrast, 
in our approach the strong suppression of heavy quarks 
is related to the steep z-dependence of the fragmentation 
function D^/ q (z). This function is fitted to data on pp 
collisions, while no fitting is done to nuclear data. 
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VIII. CONCLUSIONS AND OUTLOOK 

Nuclei suppress the large xp production of different 
species of particles: light hadrons of both small and large 
Pt, dileptons, hidden and open heavy flavor, photons, 
etc. So far, no exception is known. The source of the 
effect can be understood either as an extra Sudakov sup- 
pression caused by multiple interactions in nuclei, or as 
a nucleus-induced reduction of the survival probability 
of large rapidity gap processes. Although this effect can 
be also represented as an effective energy loss, the former 
scales in xp . This is different from a simple singlc-parton 
energy loss, which is energy independent and leads to an 
energy shift Axi = AE/E that vanishes with energy. 

In addition to this key observation, the new results of 
this paper can be presented as follows. 

• The simple formula (11) based on Glauber-Gribov 
multiple interaction theory and the AGK cutting 
rules explains the universal xp scaling observed 
in data for inclusive production of leading light 
hadrons with small pp quite well. 

• With the same input we calculated high-py hadron 
production at large xp and found a substan- 
tial suppression. This parameter-free calculation 
agrees with recent measurements performed by the 
BRAHMS collaboration at forward rapidities in 
deuteron-gold collisions at RHIC. Our simple ex- 
planation is based on just energy conservation; 
therefore, it could be implemented independently 
of the dynamics. On the other hand, it docs not 
leave much room for other mechanisms under de- 
bate, such as the CGC. We expect a similar sup- 
pression at large pt and large xp at lower energies, 
where no effect of coherence is possible. 



• The Drell-Yan process, treated like heavy photon 
bremsstrahlung in the target rest frame, is also sub- 
ject to a nuclear suppression at large xp imposed 
by energy conservation restrictions. In this case, 
we made a model-dependent calculation and relied 
on a fit to data for soft production of light hadrons. 
Within experimental uncertainties of the available 
data, we described rather well the data for DY re- 
action at large masses where an alternative expla- 
nation, nuclear shadowing, is excluded. 

• Charmonium production is different from DY only 
by a steeper fragmentation function q — > "fq, which 
we fitted to pp data, and by an additional contribu- 
tion of higher twist shadowing related to the large 
size of the charmonium. We correctly reproduced 
the magnitude of nuclear suppression at large xp, 
but the shape of the xp dependence needs to be 
improved. This problem seems to be a result of 
our model-independent, but oversimplified, fit to 
soft hadronic data. We leave this improvement for 
future work, both for charmonium production and 
the DY process. 
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